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ABSTRACT 

Saturn's magnetic field is remarkably axisymmetric. Its dipole axis is in- 
clined by less than 0.2° with respect to its rotation axis. Early evidence for 
nonaxisymmetry came from the periodicity of Saturn's kilometric radio bursts 
(SKR). Subsequently, percent variations of the SKR period were found to occur 
on timescales of years. A recent breakthrough has been the direct detection of a 
nonaxisymmetric component of the field that rotates with a period close to that of 
the SKR. Because this component's magnitude varies only weakly with distance 
from Saturn, it must be supported by currents external to the planet. These 
currents flow along field lines that connect plasma in the equatorial region of 
the outer magnetosphere (the plasma disk) to the high latitude ionosphere. The 
plasma originates from mass lost by the planet's rings and moons. It is tightly 
coupled to the magnetic field and its motion can be ascribed to the large scale 
interchange of flux tubes. Heavily loaded tubes drift outward and are replaced 
by lightly loaded ones which drift inward. This process of rotationally driven 
convection is responsible for breaking the axisymmetry of Saturn's external mag- 
netic field. The convection pattern rotates rigidly with the angular velocity of 
the plasma at its source. Its rotation provides the clock that controls the periods 
of both the SKR bursts and the nonaxisymmetric magnetic perturbations. We 
distinguish two types of currents. Those that flow in the azimuthal direction in 
both the ionosphere and the plasma disk provide a radial force that enables the 
plasma to corotate with the planet's spin. They are responsible for stretching 
magnetic field lines and thus exposing the polar ionosphere to the incoming solar 
wind. Other currents flow in the latitudinal direction in the ionosphere and the 
radial direction in the plasma disk. These act to transfer angular momentum 
from the planet's spin to the outflowing plasma. They slow down the rotation 
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rate of the ionosphere below that of the underlying atmosphere and are the rea- 
son the clock referred to above runs slow. Moreover, these currents are the source 
of nonaxisymmetric magnetic perturbations whose strength varies inversely with 
radial distance in the planet's equatorial plane. Quantitative agreement with 
the magnitude of these perturbations requires a mass loss rate of order 10 4 gs -1 , 
similar to that believed to come from Saturn's E-ring. 

1. Introduction 

Saturn's atmosphere exhibits strong (< 400 ms" 1 ) and stable (over decadal time 
scales) zonal winds. This precludes assigning a unique rotation period to its outer layers. 
Although the planet is fluid throughout, its deep interior must be in near solid body 
rotation (Liu et al. 2006). Thus the internal rotation rate might be revealed by observing 
the nonaxisymmetric components of the planet's magnetic field. 1 However, application of 
this technique to Saturn has been hampered by the extreme axisymmetry of the planet's 
magnetic field. 

Voyager observations of Saturn Kilometric Radiation (SKR) bursts coming from the 
planet's auroral regions showed a periodicity of lOh 39 min 24 ± 7 s (Desch and Kaiser 
1981). The periodicity was suspected to arise from a small non-axisymmetry of the planet's 
internally generated magnetic field. Detections of small in situ magnetic anomalies from 
Voyager and Pioneer 11 magnetometers were also reported (Giampieri and Dougherty 
2004), with periods consistent with that of the SKR. Ulysses observations of the variability 
of the SKR period, by of order 1 % on timescales of 1 year from 1994 to 1997 (Galopeau 

1 This technique applied to Jupiter has found a rotation period constant to within seconds 
over 50 years. 



-4- 



and Lecacheux 2000), challenged its interpretation as the rotation period of the planet's 
deep interior. 2 Cassini confirmed this variability, measuring an SKR period of lOh 45min 
45 ±45s on approach to Saturn in 2004 (Gurnett et al. 2005). Most recently, magnetometer 
data obtained by Cassini showed a small (~ few nT) signal with period lOh 47min 6 ±40 
s (Giampieri et al. 2006) that was stable during 14 months of observation. Giampieri et 
al. (2006) suggested that this period might be that of Saturn's interior spin. However, 
the decline of the perturbation amplitude with radius is too slow to be due to a current 
source within Saturn, so this interpretation cannot be correct. Instead, we propose that 
the non-axisymmetric component of Saturn's external magnetic field is generated by 
rotationally driven convection in the planet's magnetosphere. The convection transports 
plasma from the inner magnetosphere to the magnetopause where it joins the solar wind. 
The plan of the paper is as follows. In §2 we provide a simplified version of equations 
governing rotationally driven magnetospheric convection. We apply these equations in §3 
to estimate the nonaxisymmetric magnetic perturbations it produces. §4 is devoted to a 
discussion of the clock that controls the perturbations' rotation rate. A short summary is 
given in §5. 

1.1. Nominal parameters 

We adopt cgs units for length, mass, time, and Gaussian units for electrodynamical 
quanitites. In order to focus our discussion on Saturn, we provide numerical estimates along 
with some of the major equations. The parameters used in these evaluations are displayed 
in Table 1. 

2 A change of 1 % in spin period of the planet's core over such a short timescale is 
energetically impossible. 
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Quantity 


Value adopted 


Notes 


R 


6 x 10 9 cm 




Q 


1.6 x 1(T 4 s~ 4 




Bp 


0.2 G 




Yip 


10 13 S P)13 cm s- 1 


[1] 


a,i 


4-R Cli,4R 


[2] 








M 


10 4 M 4 g s- 1 


[3] 


n 


3 x 10- 24 g 




7 


1.4 




T 


140 T 140 K 


[4] 


V 


10 8 u 8 cm 2 s" 1 


[5] 


P 1 


3 x 10" 13 pL 12 5 g cm" 


" 3 [6] 


9 


1 x 10 3 cm s~ 2 





Table 1: The numerical quantities adopted in order of magnitude calculations throughout 
the paper. The meanings of the symbols are given in the text where first used. Key to 
numbered notes: [1] scaled between estimates for low and high latitude ionosphere (Atreya 
et al. 1974); [2] scaled to orbit of E-ring; [3] scaled to estimates for plasma production 
from Saturn's E-ring (Leisner et al. 2006); [4] scaled to atmospheric temperature at base of 
ionosphere (Atreya et al. 1974); [5] scaled to eddy diffusion coefficient at homopause (Atreya 
et al. 1974); [6] scaled between estimates for low latitude ionosphere and auroral ionosphere 
(Atreya et al. 1974). 
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2. Basic Electrodynamics 

2.1. Rotationally driven magnetospheric convection 

The topic has an extensive history, but this is not the place to review it. Instead, we 
point to a few influential papers that aided our understanding. We have done nothing more 
than to apply what we learned from reading the literature. 

Equations governing rotationally driven convection were formulated by Chen (1977) 
and by Hill et al. (1981). Many applications have been made to the outward transport 
of mass from the Io plasma torus. Pontius and Hill (1989) includes a clear discussion of 
different approaches to this problem. 

Progress in solving the equations referred to above has been slow. In retrospect, 
this is not surprising. They are nonlinear set, which probably precludes finding analytic 
solutions. Moreover, realistic applications are faced with including a continuous supply 
of plasma along with boundary conditions that simulate its loss to the solar wind at the 
magnetopause. 

Perhaps the most ambitious attempt at a realistic solution is that by Yang et al. (1994) 
who applied the Rice convection model to the Io torus. They investigated an initial value 
problem, the instability of a torus of finite width. An active source of plasma was not 
included. Long fingers were found to grow radially outward from the torus. This is not 
surprising. The initial state is analogous to that of a heavier fluid resting on top of a lighter 
one which is Rayleigh- Taylor unstable. But unlike the standard Rayleigh- Taylor instability 
which takes place for constant gravitational acceleration, the instability of the plasma torus 
is driven by centrifugal acceleration which increases linearly outward. This increase allows 
narrow fingers to run away from the more slowly developing, thicker modes that might 
otherwise subsume them. 
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In the absence of anything better, we adopt a simplified picture of steady-state 
magnetospheric convection. A tongue of plasma flows outward from a torus of neutral 
material that is escaping from Enceladus. Plasma is continuously created inside the torus 
by ionization of this material. It exits in a tongue which flows outward. Except within the 
tongue, flux tubes outside the torus drift inward due to "fringing" electric fields surrounding 
the tongue region. As they cross the torus, the inward-moving tubes are loaded with freshly 
created plasma. Then their trajectories bend around so that they join the back of the 
tongue. In this manner, plasma is continuously removed from the entire torus even though 
the tongue emanates from only a limited range of azimuth. 

In order for the convection pattern to remain steady, plasma must consistently outflow 
from the same range of azimuth in the rotating frame. Outflow occurs from the densest 
part of the torus, and so if a single tongue is to carry the outgoing material, its base must 
always be refilled fast enough so that no other longitudes in the torus can accumulate more 
plasma. Inwardly drifting tubes, which empty the rest of the torus of plasma, drift more 
slowly than those moving outwards, because the electric fields in the inward-drifing tubes 
cannot be larger than the internal tongue fields which they "fringe" . Only if the base of the 
tongue spans A<pi > n radians can the inwards drift velocity be large enough to keep the 
plasma content in the rest of the torus lower than in the tongue region. 

We have been able to bolster this description with simplified models of electrostatic 
fields, but much remains to be done before anything rigorous might emerge. A serious 
technical issue is that a smooth tongue of plasma is likely to develop narrower fingers as 
discussed above. Fortunately, the conclusions of our investigation are insensitive to this 
possibility. However, it would certainly impede a rigorous calculation of the convection 
pattern. That will have to be left for the future. 
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2.2. Notation 

We adopt spherical polar coordinates r, 9, and work in the inertial frame. Saturn's 
magnetic field is approximated as a spin aligned dipole. The ionosphere is taken to 
rotate with uniform angular velocity Q. 3 Where necessary, superscripts M and 1 are used 
to distinguish magnetospheric and ionospheric quantities, and Ip to denote the direct 
(Pedersen) component of the ionospheric current. B r and Bg are the components of the 
unperturbed magnetic field at a general field point. R is Saturn's radius and a is the orbital 
radius in the equatorial plane. B P is the magnetic field intensity on the equator in Saturn's 
ionosphere and B z is the component of the vertical magnetic field in the magnetic equator 
at a; 

*.=-©V (i) 

Height integrated current densities and electrical conductivities are indicated by J and S, 
respectively. The surface mass density in the magnetosphere is denoted by a. 



2.3. Dipole magnetic fields 

Components of a spin aligned dipole magnetic field take the form 

2Mcos9 MsinO 

where M is the dipole moment. The field magnitude is 

B=(B 2 + B 2 ) 1/2 = ^(l + 3cos 2 9) 1 / 2 . (3) 

A individual field line is labeled by either the colatitude of its footprint at r = R, denoted 
by 6*o, or by its maximum radial extent a; asm 2 9 = R. Its shape is described by: 

r sin 2 9 = R sin 2 9 or r = a sin 2 9 . (4) 



3 Except in §4 where f2 denotes the angular velocity of the deep atmosphere. 
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We are interested in field lines that connect to the planet at high latitudes. Thus we 
simplify our expressions by setting sin9 = (R/a) 1 / 2 and cos# — 1- 

2.4. Magnetospheric currents 

Consider an element of cold plasma that is nearly corotating with and slowly drifing 
away from the planet. Centrifugal balance and angular momentum conservation require 
that the height integrated current densities which pass through the element satisfy 4 

= -atfa (5) 

c 

and 

= -2and (6) 

c 

respectively. 
The ratio 



Jf _ 2a 



m (7) 



is small in the inner magnetosphere but becomes of order unity in its outer regions. 



2.5. Ionospheric direct currents 

The horizontal divergence of the magnetospheric currents, V2d • J M , flows along 
magnetic field lines and closes in the ionosphere. This determines the components of the 
direct (Pedersen) current in the ionosphere. They read 

T i p 1 da Jf can 2 R /a\H/2 



Rd6 n 2 Bp \RJ 



4 We neglect the planet's gravity in equation (5) 
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and 

jip _ 

" Rsm6 2 Bp \R 



ji P = a Jg = cvQR_ / ayi/2 d (g) 



Their ratio is given by 

(10) 



The ionospheric Hall current, J 7/l , is not fixed in this manner. Because its horizontal 
divergence vanishes, 5 its determination requires knowledge of the ionospheric electric field. 



2.6. Ionospheric electric field 

The height-integrated current density, J 7 , and electric field, E 7 , are related by 6 

J 7 — J*p + J Ih — E P F 7 + E/j(b x E 7 ) , (11) 

where F 7 = E 7 + (Q x R) x B 7 /c is the Lorentz force, S p and are the height-integrated 
Pederson and Hall conductivities, and b = B 7 /|B 7 |. Taking the divergence of equation (11) 
yields 

V M • J 7 = V 2d ■ 3 Ip = S P V 2 , • E 7 . (12) 

V 2 d ■ J Ih = because V 2 d x E 7 = 0; E 7 is a potential field. A full solution for E 7 would 
involve setting E = — V$ 7 and then solving Poisson's equation V 2 $ 7 = — S p V2d • J 7p - As 
detailed in §2.1, for realistic conditions, this has proven to be a difficult task. Fortunately, 
a simpler procedure suffices for the purposes of the current investigation. 

We set 

jip jip 
Fe = ^- and ^ = -J- , (13) 

5 Provided the Hall conductivity is independent of position as we assume it to be. 
6 By setting sin# = (R/a) 1 ^ 2 and cos#o = 1, we neglect both B\ and the effects of the 
parallel conductivity. 



- 11 - 



with 



2ORsin0 „ . , 

Fg = E e Bp. (14) 



This procedure does a good job evaluating E 7 in the portion of the ionosphere that is 
magnetically connected to the outgoing tongue of magnetospheric plasma. However, it does 
not permit a determination of the fields that fringe this region. These control the inward 
flow of depleted plasma tubes. 



2.7. ExB drift 



Just above the ionosphere, the plasma drifts at velocity 



V = C( • (15) 



Projecting down to the magnetosphere, we obtain 

a 26 cEs /a\V2 c 2 aVt 2 fa 



and 



6q RBp \RJ YjpB'jj \R 



An cF e ( a \ 1/2 _ c 2 a ( a \ 6 a 



(16) 



, , (17) 

n 2QRB P \RJ 2H p B 2 p \r) a' K J 

Here Afl is the angular velocity at which the plasma in the tongue slips relative to the 
rotation at the top of the ionosphere. 

For future reference, we note that 

' \ 2 

2nAtt. (18) 
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2.8. Coupling to rate of mass loss 

Suppose the tongue covers A<f>(a) in azimuth where A0 is a function of a to be 
determined later. Then 

M = A<paaa. (19) 

Substituting for a using equation (16), we arrive at 

• A<j)c 2 a 2 n 2 R 2 /a\8 . . 

~ ^ y • < 2o > 



from which we obtain 

°" ~ dlR ^ A0 
Next we replace cr in a/ a which yields 

, 1/2 

a Vtc I M \ /a\ 2 



a R\B P \ \ S p A0 / Vi? 



(21) 



(22) 



2.9. Steady-state scalings 

Flux freezing implies that a/B z is independent of a in a sourceless, steady-state flow. 
Consequently, A<p at a is related to its initial value at Oj by 

A0=(^) 2 A0,. (23) 



Thus 



cfiaj \ A0j / \ a 



* V2 

i?fic I M \ /a\ 4 



(24) 



o • (25) 



Somewhat arbitrarily, we adopt a , the value of a where a = Qa, or equivalently where 
AQ = — 0/2, as the outer radius of the region in which partial corotation applies; 

p j = ^\Bp\_^ 1/3 ^A^j 1/6 ^ 21 ^ a|^A^S^ ^ 1/6 (2g) 
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The timescale a/ a for outward plasma transport decreases with distance from the 
planet. Plasma near the source at Oj doubles its radial distance in a time ~ 10 days. 
The constancy of the observed magnetic period on far longer timescales implies that the 
convection pattern remains steady for many dynamical timescales of the source region. 



3. Magnetic Perturbations 

3.1. Opening of the polar cusp 

The azimuthal current J^f creates a radial component of magnetic field just above and 
below the tongue of outgoing plasma; 

2nJ M 

AB r = ^sgn(^) . (27) 

c 

Applying equations (24) and (26), we find that these field lines bulge outward relative to 
vertical by an angle 

B, Bl c' \aj \ M, J 1 ' 

The stretching of these field lines is a consequence of the centripetal acceleration they 
impart to the plasma. 

If a attains a substantial value at a Q , the exposure of Saturn's polar cap to the incoming 
solar wind will be enhanced along the range of longitudes subtended by the outer parts of 
the tongue. 7 This could account for the preferential emission of the SKR within a narrow 
range of longitudes. 

7 The field lines bulge only above and below the outgoing plasma. Because A</> oc cr 2 , 
A0 O m &y be quite small, even if A0» is of order n. 
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3.2. Magnetic anomalies measurable by spacecraft 

The radial current in the tongue of outgoing material at a is given by 

I? = aA0jf = -^L . (29) 

Thus 

rf/ a M , 6cQM , . , 

— ^da = — da 30 

da aB z 

is the current that flows along the strip of field lines connecting the ionosphere (in each 
hemisphere) to the tongue between a and a + da. It is notable that iff depends on the rate 
of mass loss but not on the pattern of outflow. 

To estimate the magnetic perturbation produced by this current, we approximate it as 
flowing along a thin wire of infinite extent which passes the observation point at distance s. 
Then 

d5B 2 dljf 12fiM/a\2 

« — ^ — . (31) 

da cs da RB P s \RJ 

We focus on perturbations in the equatorial region of the inner magnetosphere, since 
this is where the Cassini magnetometer has found evidence for non-axisymmetric field 
components. For an observer in the equatorial plane at radius zu, the minimum distance s 
between the observer and a "wire" of dipole field-aligned current which passes through the 
equator at a ^> zu is s ~ zu. Substituting for s in Eq. 31 and integrating to a , we find 



5B 



f^F (I)' w ( E ^ A ^) V2 ~ 12 x 10 " 5 G (^>4S P>13 ) 1/2 • (32) 



The above estimate for the magnitude of the magnetic perturbations ignores the 
contribution from the return current that comes from a > a Q and closes the circuit. 
Unfortunately, not much can be deduced about the geometry of the return current since 
the field lines it flows along are likely to be strongly perturbed by the solar wind. However, 
it is possible that it could act to reduce the perturbation magnitude, perhaps by a factor 2. 
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We have also ignored the contribution from the currents that flow in the opposite (N versus 
S) hemisphere. In the case of perfect N-S symmetry, magnetic perturbations from the two 
hemispheres would exactly cancel in the middle of the tongue but not elsewhere. In reality, 
at most times 8 this cancelation is likely to be small, since the ionospheric conductivity will 
differ in the northern and southern polar ionospheres. Near solstice, almost all the current 
will flow through the summer hemisphere. At these times the estimate for SB should be 
doubled. 

As illustrated in Figure 1, the field-aligned (Birkeland) currents supplying J^ 1 to the 
tongue occur in closely spaced, oppositely directed pairs. The magnetic perturbations 
produced by these "line dipoles" are small at distances large compared to their spacing. 
Thus, although they are typically smaller, the field-aligned currents feeding dominate 
the magnetic perturbations present in the magnetosphere. 

The fragmentation of the plasma tongue into narrower structures would have a minor 
effect on the magnetic perturbations produced by the currents supplying j£* , but would 
further weaken those produced by the currents supplying J™ . 

Analysis of the variation of the measured components of 5B with position in the 
magnetosphere would test our model. 

4. The Clock 

As described in §1, a somewhat slow and imperfect clock controls the quasi-periodic 
behavior of SKR bursts. In our scenario the same clock controls the magnetic anomalies. 
Both Voyager- and Cassini-eia. measurements of SKR and magnetic periods are consistent 



The exception being near times of equinox. 
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phi-currents close in 




Fig. 1. — The currents associated with the plasma tongue in Saturn's magnetosphere: az- 
imuthal currents occur in close opposite pairs, and so their contribution to magnetic pertur- 
bations is small. 

with this picture. 

4.1. Where is the clock located? 

The clock is located in the inner magnetosphere and beats at the period of rotation 
of the asymmetry associated with the inner portion of the plasma outflow. 9 As explained 
below, the clock's period propagates throughout the magnetosphere provided the differential 

9 Our working hypothesis is that the E-ring centered at a ~ 4R around the orbit of 
Enceladus is the dominant plasma source. 
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rotation of the latter remains time invariant. 

Suppose that a source of material located at radius feeds a plasma torus rotating at 
The torus is unstable and sends out a "tongue" of plasma centered on a fixed azimuth 
4> c {a>i) in the rotating frame. Plasma at the center of the tongue moves outwards at radial 
speed a c (a), and orbits at angular speed </> c (a), where 4> c (a) < Jlj. In steady-state, the shape 
of the tongue's centerline would be determined by 

fda'^-. (33) 

Viewed from a nonrotating frame, the tongue is a steady structure rotating at pattern 
speed fli. At radius a and time t, the apparent azimuth of the tongue's centerline is given 
by 

0obs = <pc(a) + ^it, (34) 
and so o b s = regardless of the run of differential rotation across the magnetosphere. 



4.2. Why is the clock slow? 

SKR and magnetic periods determined by Cassini are longer than those associated 
with the motion of any atmospheric features. The most plausible explanation is that the 
ionosphere rotates more slowly than the atmosphere below it because magnetic torques are 
transferring angular momentum from it to the plasma tongue. 10 

Since the magnetic torque increases sharply with increasing latitude, this is also a 
plausible explanation for the observed decline of the magnetosphere's angular velocity with 

10 Huang and Hill. (1989) concluded that the rotation of Jupiter's ionosphere is slowed in 
this manner. 
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increasing distance from Saturn. We shall show that this scenario is qualitatively reasonable, 
in contrast to models based on the slippage of the rotation of the magnetospheric plasma 
relative to the rotation of the ionosphere. Most of the way out, the plasma maintains good 
corotation with the part of the ionosphere to which it is connected, but the ionosphere is 
subcorotating with respect to the underlying atmosphere. The observed clock frequency is 
the rotation rate of the ionosphere where it connects to the inner part of the plasma tongue. 



4-2.1. Steady-state rotation of the ionosphere 

We analyze a simple model for the steady-state rotation of the ionosphere and 
underlying atmosphere. It assumes axial symmetry and considers only vertical transport 
of angular momentum. Deep atmospheric layers are taken to rotate rigidly with angular 
velocity Q. We work within the approximation of an isothermal atmosphere with sound 
speed c s , scale height H = c^/jg, buoyancy frequency N 2 = (7 — l)g/(jH), and eddy 
diffusivity v. The ionosphere is taken to be a single layer rotating at the angular velocity 
Vl M of the part of the plasma tongue to which it is magnetically connected. 

We modify equation (6) to allow for the nonuniform rotation rate, f2 M , of the 
magnetospheric plasma. The torque per unit a applied to the tongue of outgoing plasma 
reads 

dT M d 

= M-— (n M a 2 ) . (35) 
da da 

Thus the magnetic torque per unit 6q on the northern ionosphere is given by 



dT^ MR 2 d fVL 



(36) 



d0Q 2 d0Q \ 9q 

In steady state, the torque must be constant with depth below the ionosphere. Provided 
the torque is not too large, a stable, steady-state, shear flow is established in which the 



-19- 



viscous torque given by 

dT v 3 dn A 

— = -2nR9 oP v—, (37) 

carries angular momentum up from the atmosphere to the ionosphere. Equating the viscous 
torque to the magnetic torque yields 

dfl A _ M d /fi M \ 

~dz~ ~ AntVpvOl df [if ) ' (38 ^ 

The shear flow is stable where the Richardson criterion is satisfied, that is where equation 
(38) predicts 

Mo^<N. (39) 

In the stable regime, 

q m MH d (Q M \ MVLH 

Here we have set Vt M = VL in the final step. This is a good approximation since Vl — Vl M <^ VL 
in the stable regime. 

The boundary of the stable regime occurs where the stability criterion is violated just 
below the ionosphere. We use the symbol O to denote the value of 9 at this boundary. At 

e , 

m d fn M \ 

~ -1 . (41) 



iTiRp^NQl d9 V 9^ 

An approximate solution for O , and the corresponding a crit /R = 0q 2 , follows from setting 
Q M = Q: In this manner we arrive at 

2/7 / i \ 2/7 



^crit 



The unstable layer penetrates deeper into the atmosphere poleward of 0o- Because the 
angular velocity gradient in the stable layer is inversely proportional to density, and the 
density increases exponentially with depth, the thickness of the unstable layer increases 
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logarithmically with decreasing 9 . An estimate for ft — ft M is obtained by multiplying the 
critical angular velocity gradient from equation (39) by the thickness of the unstable layer. 
The following expression provides a good fit to ft — ft M for all values of 9 or a/R. We 
express it in terms of the latter for ease of comparison with data on the rotation of plasma 
in Saturn's magnetosphere. 



ft M NH / a\ 1/2 

1 ~ — In 

ft SIR \rJ 



1 | 7/2 a dVt M 



a C rit 7 V ft 2ft da 



(43) 



Figure 2 displays the run of = ft M a vs a/R for our nominal parameters. The 
rotation curve is in reasonable agreement with Voyager measurements (Richardson 1986). 
At a/R = 4, corresponding to the orbit of Enceladus and the brightest portion of the 
E-ring, 1 — ft M /ft ~ 0.5%. By comparison, equation (18) predicts the much smaller value, 
3 x 10~ 5 , for the slippage of the rotation rate in the plasma tongue at a = AR relative to 
that of the ionosphere. 



4.3. Why Is the Clock Imperfect? 

The ionospheric rotation rate responds to changes in a variety of parameters, including 
the atmosphere's eddy diffusion coefficient, the ionosphere's height-integrated conductivity, 
and the rate of mass loss from the magnetosphere. All of these are likely to have a seasonal 
dependence. Solar weather, including short term variations in the solar wind ram pressure 
and longer term variations over the solar activity cycle also affect the conductivity in the 
auroral ionosphere. 
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Fig. 2. — Orbital velocity of equatorial plasma as a function of distance from Saturn. Rigid 
corotation with Saturn's interior is plotted as the dashed line; the magnetosphere lags the 
planetary interior due to slowing of the ionosphere by magnetic torques. 

5. Summary 

Rotationally driven convection of magnetospheric plasma breaks the axisymmetry of 
Saturn's external magnetic field. Field aligned currents transfer angular momentum from 
the planet to a tongue of outflowing plasma. This transfer slows the rate of rotation of 
the ionosphere relative to that of the underlying atmosphere. The currents are the source 
for the non-axisymmetric components of the field. The common rotation rates of these 
components and Saturn's kilometric radio (SKR) bursts is that of the plasma near the orbit 
of Enceladus, and by extension the rotation rate in the ionosphere to which this plasma is 
coupled. This rate tells us nothing about the rotation rate of Saturn's deep interior. Of 
that we remain ignorant. 
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Magnetic perturbations with magnitudes similar to those observed by Cassini are 
produced for M m 10 4 g s™ 1 , a value similar to estimates for the rate of production of 
plasma from Saturn's E-ring. 

Enhancement of the SKR occurs in a narrow range of longitudes where the tip of the 
outgoing plasma stream connects to the auroral ionosphere via field lines that are bowed 
outwards by currents that supply the plasma's centripetal acceleration. 
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